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SUMMARY 


The computer code " AGDISP" ( Agricultural DIS persal) has been developed to 
predict the deposition of material released from fixed and rotary wing 
aircraft in a single-pass, computationally efficient manner. The formulation 
of the code is novel in that the mean particle trajectory and the variance 
about the mean resulting from turbulent fluid fluctuations are simultaneously 
predicted. The code presently includes the capability of assessing the 
influence of neutral atmospheric conditions, inviscid wake vortices, particle 
evaporation, plant canopy and terrain on the deposition pattern. 

The AGDISP code is appropriate for predicting the motion of material 
released for times over which aircraft-unique wake and propulsion system 
effects are still influencing the motion of the released material, and 
sufficient diffusion has not resulted in a distinct spray cloud. For later 
times Gaussian plume modeling or other methodologies may become more 
efficient. AGDISP output is therefore configured in a format which may be 
used as input to these later-time codes. The AGDISP code has an additional 
flexibility in allowing the user to specify a flow field independent of the 
flow model options available in the code. This feature gives AGDISP the 
capability of investigating particle dispersion in specific situations of 
interest. 
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1 . INTRODUCTION 


The development of "AGDISP ,f was motivated by a desire to determine how 
aircraft-unique wake and propulsion system characteristics affect the ground 
deposition pattern of aerially released material. Development of aircraft 
flow field models during NASA’s vortex wake hazard program provided 
sufficiently detailed, yet simple, models (at least for fixed wing aircraft) 
which could be utilized in this application. What was needed was the 
development of an efficient algorithm permitting the prediction of particle 
dynamics including the effects of turbulent fluid fluctuations. Development 
of this algorithm as well as a description of the flow field models programmed 
into the AGDISP code form the subject of this report. 

Historically computers have been used since 1953 to compute the motion of 
material released from aircraft, with the first published study undertaken by 
Reed.^ ' While Reed grossly oversimplified the aircraft wake flow field by 
neglecting roll up of the vortex sheet trailed from the wing and omitting 
propeller swirling wash, he did acknowledge the importance of the vortex wake 
in establishing particle trajectories. The model developed by him integrated 
the equations of motion governing the dynamics of single particles in a vortex 
wake flow field modeled as two counter-rotating irrotational line vortices 
with separation distance equal to the span of the wing. Image vortices were 
used to simulate an inviscid ground plane. Resultant trajectories were 
determined for several particle sizes released from several positions along 
the trailing edge of the wing; then, assuming a particle size spectrum, an 
expression for the deposition on the ground was developed. 

The strength of the Reed model rests with its simplicity. It has been the 
starting point for several more recent formulations of the problem. Trayford 
and Welch' ' added the effects of propeller wash, cross-wind and evaporation 
to the Reed model. Their calculations showed that both the propeller wash and 
cross-wind could have a significant effect on the deposition pattern. 
Unfortunately, we believe that they did not include the effect of evaporation 
correctly in the droplet dynamic equations. When mass is lost from a droplet 
in relative motion with a gas, both the droplet and gas must feel opposite 
forces of magnitude V re jdM p /dt where V rel is the relative velocity and 

dMp/dt is the rate of mass lost. The Trayford and Welch model did not 

include this force contribution to the particle equation of motion. 
Additionally, as in the Reed model, the interaction of the vortices with the 
ground plane was assumed to be Inviscid. 

More recently, Bragg^ attempted to improve the flow field in the Reed 
model by including the effect of bound circulation on the wing as well as 
allowing for some three-dimensionality in the vortex wake. Since this model 
was used to help guide definition of subscale deposition experiments at the 
Langley Vortex Research Facility, the effect of inviscid tunnel walls was also 
included. Recently, Loats Associates, Inc.' ' reviewed additional deposition 
models used to predict aerial spray drift. These additional models, in 
general, utilize Gaussian plume modeling and thus do not give accurate 
predictions when strong aircraft flow fields are influencing particle 
behavior. 
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In a series of reports, Wickens^ 5,6 ^ discussed the aerodynamics of wakes 
from fixed wing aircraft and their implication with regard to aerial 
application. Wickens acknowledged that the vortex flow field and its 
interaction with the atmosphere and the surface is far more complicated than 
that which could be modeled by two inviscid point vortices. Recently, Bilanin 
et al. ' ', using NASA Langley’s WAKE computer code, demonstrated that the 
viscous interaction of a vortex pair with the ground results in the shedding 
of secondary vorticity from the surface. This secondary vorticity completely 
alters the trajectories of the vortices and results in the vortices actually 
moving away from the surface. Jordan et al.^®', using vortex trajectory data 
measured in the Langley Vortex Research Facility, compared measurements with 
computations of the viscous interaction of the vortex pair with the ground and 
found favorable agreement. Further support that complicated interactions may 
be expected in a canopy with and without cross shear is given in computations 
presented by Morris^. Here the descent of a vortex pair into a forest 
canopy was simulated, and the computation showed that the canopy would halt 
the lateral motion of the pair. Computations of vortex pairs in proximity to 
the ground with cross-wind shear produce a tipping of the pair, an observed 
phenomenon. 

The computer code WAKE, developed by Teske^^^^^, is capable of computing 
the turbulent evolution of an aircraft wake flow field interacting with the 
atmosphere and the ground. In principle this code can provide the flow field 
upon which a dispersal code is based. The computer resource necessary to 
generate these flow field results, however, can be quite significant. 

It would seem, therefore, that the largest obstacle to predicting the 
dispersal of material from agricultural aircraft is the determination of the 
flow field acting on the particles. In principle, however, the technologies 
required to define these turbulent flow fields are readily available. When 
the particle concentrations are sufficiently dilute, momentum transfer from 
the particle to the fluid may be neglected (the dynamics of the fluid are then 
uncoupled from the particle). What remains to be determined is whether the 
particle dispersion can be adequately described by an Eulerian or Lagrangian 
formulation within reasonable computational limits. 

This report is organized as follows. In Section 2 the equations governing 
the motion of aerially released particles are developed, including a 
description of the evaporation model used. In Section 3 the flow field 
options available in the AGDISP code are reviewed. Code limitations and 
applicability are included in Section 4 and case studies form Section 5. The 
equivalent Gaussian distribution is discussed in Section 6. Lastly in 
Section 7, conclusions and recommendations are offered. 
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2. PARTICLE DYNAMIC MODELING 


Particle Description 


The motion of particles may be computed from either an Eulerian or 
Lagrangian formulation. After analysis of computer time requirements, it was 
judged that the Lagrangian formulation would be far more economical, since 
many parametric engineering computations might ultimately result. The 
Eulerian or field description of a spray cloud would involve the numerical 
solution of partial differential equations which presently requires a large 
computer resource. On the other hand, a Lagrangian formulation, tracking the 
motion of each particle separately, requires the straightforward numerical 
solution of ordinary differential equations. For these reasons a Lagrangian 
formulation was chosen as the particle description in the AGDISP code. 


Equations Governing Particle Dynamics 

The significant forces acting on aerially released material include: 

a) Weight 

b) Aerodynamic drag 

c) Force resulting from evaporation 

As shown in Figure 1 , these forces result in an acceleration of a particle of 
mass, Mp , according to: 



«i + V -- 
p 


(U 1 + 


V J - vj dM 
I i p . 

M dt *i 

P 


viscous 

drag 


acceleration resulting 
from evaporation 


gravity 


( 1 ) 


Here, and are the ensemble averaged i t ^ 1 component of particle 

position, velocity and fluid velocity, respectively, while , v^ and u^ 

f. 1. 

are the fluctuating i Ln component of particle position, velocity and fluid 
velocity, respectively. The particle velocity is related to the particle 
position through: 


d(X i + x ± ) 
dt 


V + v 
i i 


( 2 ) 
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g 


X. , V. , U. - i ensemble averaged component of particle position, 

velocity and fluid velocity, respectively. 

t 

x. , v. , u. - i fluctuating component of particle position, 

velocity and fluid velocity, respectively. 


Figure 1. Particle notation convention 



The drag force on the particle is computed from: 
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(U. + u.- V - v.) 
1 x i l 


( 3 ) 


where C D is the Reynolds number dependent drag coefficient of the assumed 
spherical particle of diameter Dp and density Pp , and P a ^ r is the 
density of air. 

Currently, we are using a semi-empirical relationship for suggested 
by Langmuir and Blodgett^^: 


C D = H (i + 0.197Re 0,63 + 2.6 x 10 4 Re 1,38 ) 


(4) 


where Re is the Reynolds number defined as: 


Re 


Yrel 

V air 


(5) 


Equation (4) is valid for spherical particles in the range 
0 < Re < 50000 . Substituting for C D in Eq. (3) yields: 


TT -T- (u i + “i - v i ‘ V 
P P 


( 6 ) 


where is the particle relaxation time computed from: 


t = 


d 2 p 
P P 


p 18 v . p 
r axr ai 


r ( l + °‘ 


197Re°* 63 + 2.6 x 10 _4 Re 


1.38 


(7) 
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The particle relaxation time given in Eq. (7) has physical significance with 
regard to dispersion in that it is a measure of the time required for a 
particle to catch up to the local fluid velocity. Specifically, if a particle 
were released at rest in a stream of velocity! U|> the time required for the 
particle to reach 0 * 6 3 1 U | is . 

If Eq. (1) and Eq. (2) are ensemble averaged with the drag model from Eq. 
(6) and fluctuations in the evaporation rate are ignored, there results: 


d 2 X. 

l 

dt 


i i dM 

• <U i- V [- + M-dI £ ] + *i 


( 8 ) 


dX l 

dt 


= V, 


(9) 


Subtracting Eq. (8) from Eq. (1) and Eq. (9) from Eq. (2) results in equations 
for the fluctuations: 





V 



. dM 

i E 1 

M dt J 
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( 10 ) 


dx. 


dt 


= v. 


(ID 


which may be premultiplied by x^ and v^ , ensemble averaged and manipulated 
to yield: 


a < x 1 x 1 > * 2<x i v i > 


(12) 


1 l dM 

<x v > = (<x u > - <x v >) [— + -i- — E ] + <v.v.> 

dt ii ii ii L x Mdt J 11 

P P 


(13) 


j , dM 

- 7 - <v.v.> = 2(<u v > - <v v.>) [ — + — E 1 

dt ii ii i i ' l t M dt J 


(14) 


P P 



where no summation is implied* The symbol < > denotes the ensemble average 

of the quantity. 

Equations (12), (13) and (14) are nine equations for the fifteen 

variables ^ x i v i^ » ^ v i v i^ * ^ x i x i^ * ^ x i u i^ and <u^v^> • Closure is 

completed by specifying approximate relationships for the quantities <x^u i > 

and < ^ u i v i" > » the correlations of the particle position and velocity 
fluctuations, respectively, with the mean velocity fluctuation. 

These relationships are developed by assuming the fluctuating fluid 
velocity is given by: 



iu>t 

(cj)e du) 


and solving 


dVf (u - v ) 


dt T 

P 



(13) 


(16) 


for v^ and x^ , then forming the products x ± u ± anc ^ u i v i anc * ensemble 
averaging. With zero initial conditions, Eqs. (13) and (16) yield: 
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(17) 


(18) 


Upon multiplying Eq. (17) and Eq. (18) by and ensemble averaging, we 

obtain: 


0 -t/x 

x e 1 
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and 


<u iV 


-t/x 

$ (oj) 1 - e ^ cos u)t 
u 


-t/T 

- e ^ x a) sin 0)t 
P 


>doj 


( 20 ) 


where $ u (w) is the spectral density function for transverse velocity 
fluctuations. The spectral density function will be selected so that Eqs. 
(19) and (20) can be integrated and simple closure achieved* 

Von Karman and Howarth^~^ showed for isotropic turbulence that 
fluctuations normal to the mean flow may be expressed as: 


$ 

u 


(a)) = A— A 

V } 3tt U 


2 1 + 3(ujA/U) 2 

9 . 9 9 

[l + (wA/U) ] 


( 21 ) 


where A is the integral scale of the turbulence, is the mean square 

turbulence level and U is the free stream velocity. If U were now 
interpreted as the relative mean velocity between the aerially released 
material and the mean fluid velocity, this form of $ u (m) permits an 

analytical evaluation of <x jL u i > and <u i v i > . Implicit in this substitution 

is the assumption that the aerially released particle is acted upon by a 
locally isotropic turbulent field (this assumption is compatible with the 
simple flow field models developed in the next section). Integration of Eqs. 
(19) and (20) with the assumed form of $ u (w) given in Eq. (21) yields: 


2 

<x.u> = 
l i 3 


-T 



( 22 ) 


2 

<u . v . > = -3- K 

li 3 


(23) 


with 


I 


8 



(24) 
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where t t = A/|V re ^| and is the time of travel of the material through the 

turbulent eddy of scale A . When the particle adjustment time T « t t , 
the particle tracks the fluid motion and: 


<x . u > 
1 i 



(25) 


<u v.> 
i l 


2 

_ 9 _ 

3 


(26) 


Assuming that ensemble averaged fluid velocities are known and fluid 
fluctuating variances are prescribed, the formulation of the particle 
dispersal dynamics problem is complete once initial conditions are specified 
at t = 0 ♦ 


Typically: X i = X io 

V i = V io 

<Xi x i> = <x iV .> = 0 
<v i v i > = <v i v i > 0 


(27) 


A non-zero specification of <v^v ± > allows for variance in the initial 

particle/droplet velocities at the point of release. The above equations are 
programmed in the AGDISP computer code. 

To verify that the formulation of the equations governing particle 
displacement variance is correct, we assume that the particle has zero mass 
and as a result passively follows the fluid fluctuations. Then x p -► 0 and 

we may determine <x^x i > for time t > x . Neglecting evaporation, Eqs. 
(14) and (26) yield: 


<v^v^> = <u 


iV 


2 

3 


(28) 


9 



Substituting this result in Eq. (13), using Eq. (25), yields: 


d__ 

dt 


<x v > = — 
i i x 

P 



<x . V > 

i i 



( 29 ) 


which for large t 


may be solved to give: 


<x i v i > 



(30) 


Therefore, from Eq. 


(12) we obtain: 


<x^x^> = 


1 q 2 t 


(31) 


This result may be compared with the result obtained using 
closure turbulent modeling for a passive tracer^^ in the t > t t 


second-order 
limit : 


<x i x i> = f 


qAt 


(32) 


Hence, when there is no relative ensemble average velocity between the 
particle and the fluid, t t must reduce to 8A/3q ♦ Therefore, the turbulent 

time scale t t must be redefined to be: 


T 

T 


A 

( l V rell + I "> 


(33) 
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Evaporation Model 


The force on a particle as a consequence of evaporation results from the 
evaporating fluid leaving the droplet at a velocity + v^ which is not in 

general equal to the local fluid velocity U<^ + u^ ; hence, a momentum 

exchange occurs. For the present a simple evaporation model v J is assumed: 


M 

P 


dM 

_£ = 
dt 




where is the life of the droplet and is given by: 


T 


e 



3A0 


t 


= 0 


( 34 ) 


( 35 ) 


where AO is the wet bulb depression and (3 


is defined as: 


6 = 84.76 [ 1.0 + 0.27 Re 1/2 ] 


o 

in units of m /sec - °C . 


( 36 ) 


During computation a particle is permitted to evaporate to a specified 
minimum particle diameter, at which point the particle diameter is held 
constant for the remainder of the computation. The minimum particle diameter 
is determined by estimating the point at which sufficient solvent has 
evaporated so that evaporation is inhibited by the high concentration of 
solute. 

Recently, Dennison and Wedding^*^ have conducted evaporation experiments 
at Colorado State University. A number of pesticides were held in a water 
solvent in a wind tunnel at terminal velocity using a very thin wire which 
spanned the tunnel. The current evaporation model was checked against this 
data. Predictions of droplet diameters versus elapsed time for the initial 
droplet size of 110 and 400 ym are shown in Figures 2 and 3, respectively. 
Solution 12 is pure water. The test data and simulation were undertaken at a 
nominal air temperature and relative humidity of 20 deg C and 20%, respec- 
tively. This corresponds to a wet bulb depression of about 10.7 deg C. As 
can be seen the predictions of particle diameter are in favorable agreement 
with test data. 
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Solutions 12 to 15 
Temperature = 20°C 
Humidity = 20% relative 
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Figure 2. Prediction and measurement of droplet diameter variation with time as a consequence of 


Solutions I2 t l3,l4 
Temperature = 20°C 
Humidity s 20% relative 
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Figure 3. Prediction and measurement of droplet diameter variation with time as a consequence of 
evaporation 


3. FLOW FIELD MODELING 


The AGDISP code has been configured to accept user-specified mean velocity and 
turbulence flow fields in neutral atmospheres. Input file format requirements are 
detailed in the AGDISP user manual ^ 1 \ For users not wishing to supply their own flow 
field, the following models have been incorporated into the AGDISP code. 

1) Aircraft vortex wake including wake roll up 

2) Helicopter flow field 

3) Propeller swirl 

4) Simple terrain 

5) Cross-wind 

6) Plant canopy 

7) Superequilibrium turbulence 

The remainder of this chapter discusses these models. 


Aircraft Vortex Wake Including Wake Roll Up 

It is well known that fixed wing aircraft trail vortex wakes as a 
consequence of the aerodynamic lift generated by the wing. This wake flow 
field can significantly influence the dynamics of material released into it. 
For this reason it is important to have available in the AGDISP code a flow 
field model which reasonably represents the velocity field in the vortex wake 
of an aircraft releasing agricultural material. Fortunately, recent research 
sponsored primarily by NASA has resulted in an understanding of the 
aerodynamics of vortex wakes (the interested reader may refer to Ref. 17 for a 
more complete discussion of this subject). The following discussion 
highlights the necessary details required to understand the flow field 
description implemented in the AGDISP code. 

When an aircraft flies at constant altitude, the aerodynamic lift, L , 
generated by the lifting surfaces of the aircraft equals the aircraft weight, 
W . Since the majority of the lift is carried by the wings and not by the 
tail, it is common practice to neglect the contribution to the wake of the 
tail aerodynamics. The wing lift distribution generates one or more pairs of 
swirling masses of air downstream of the aircraft (see Figure 4). These 
swirling masses are known as vortices. The strength of a vortex is denoted by 
the circulation r , and is related to the lift by: 


L = 


p . u br 

air o“ 


(37) 
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Vorticity is trailed forming a sheet 
directly behind the aircraft 
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Figure 4. A schematic of the roll up of a vortex sheet trailed from a simply loaded wing 




















where U^, and b are the flight speed of the aircraft and separation 
between the vortices, respectively. Unfortunately, the vortex wake 
immediately behind the aircraft (as shown in Figure 4) is generally much more 
complicated since the wake is in fact a thin sheet of vorticity trailed from 
the wing. This sheet rolls up downstream into discrete vortex pairs. The wake 
shown in Figure 4 is for a simply loaded wing; that is, a wing with a load 
distribution which trails a sheet of vorticity that will roll up into one pair 
of counter-rotating vortices. The wake flow field resulting from a lifting 
wing will become more complicated when this simple load distribution is 
altered by a deployed flap, as shown in Figure 5. In this case, in the far 
wake, two vortex pairs eventually roll up and descend under their mutual 
influence. 


In the AGDISP code the far field swirling velocities in the vortices are 
estimated using an approximate methodology first suggested by Betz.'^®' This 
methodology relates the swirling velocity distribution in the vortices to the 
details of the wing spanwise load distribution by assuming that angular 
momentum is approximately conserved. Figure 6 shows the swirling velocity 
distribution which results after roll up is complete behind wings that have 
spanwise load distributions that are triangular, elliptical and rectangular. 
The details of how these swirling velocities are computed from the spanwise 
load distribution is documented in Ref. 19. 


The above discussion presumes that the aircraft is at a sufficiently high 
altitude such that wake interaction with the ground may be neglected. Under 
most conditions of interest with regard to aerial application, however, this 
is not a good approximation and the influence of the ground must be 
included. To first order the most significant effect of the ground is to 
force the vertical component of fluid velocity to be zero at the surface. 

Out of_ground influence, the vortex pair descends downward with a velocity 
of r/2iTb while maintaining a constant separation distance between the 
vortices. The most important near ground effect is to force vortex motion 
along the surface as shown schematically in Figure 7. This effect of the 
ground plane on vortex motion is included in the AGDISP code. 

Since all agricultural material must for practical reasons be released in 
the immediate vicinity of the aircraft, the wake flow field model must be 
extended, if sufficient predictive accuracy demands it, to properly account 
for the unrolled portion of the aircraft wake. The details of treating sheet 
roll up into discrete vortices is documented in Ref. 19. Here we will 
highlight the general features of the Betz roll up model included in the 
AGDISP code. 


Figure 5 illustrates the roll up of a tip and interior flap vortex. The 
vortical material trailed from the trailing edge of the wing makes up a sheet 
which is convected about the local concentrations of vorticity. These local 
concentrations are located spanwise at positions where the slope of the wing 
spanwise load distribution is a local maximum (in practice these positions 
occur where edges of flaps or spoilers are located). Downstream of the wing, 
the AGDISP code rolls up the trailed vorticity in a manner analogous to the 
roll up of a window shade. The strength of the vortex increases as additional 
vorticity is brought into the rolled up vortex at velocity 
U t (t) = Tj./2ttRj. ( t ) for a tip vortex. The circulation increases according to: 
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Elliptical 
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Elliptical 
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Figure 6. Swirl velocity distribution of a fully rolled up wake as a function 
of spanwise lift distribution and as measured from the center of the 
resulting vortex 
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Figure 7. Trajectories of a vortex pair approaching the ground 
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dr 


t 


dt 



d£ 

dy 


( 38 ) 


where dT/dy is the strength of the vortex sheet rolling into the vortex 
and R t (t) is the Betz vortex radius computed as described in Ref. 19. 

While the vortices roll up, they are also in motion relative to each 
other. It may be shown that for two-dimensional vortices the motion of the 
centroid of a vortex is given by: 


-► 

dr 

c 

dt 


(V + U 3 k)t 


dydz 


(39) 


where the centroid is defined as: 


-► 
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(yj + zkK 

r 


dydz 


(40) 


and C is the vorticity. The velocity components Uo and U-j are the 
velocities resulting from all additional concentrations of vorticity excluding 
the centroid under consideration. During the evolution of the wake flow 
field, the AGDISP code tracks the position of the discrete concentrations of 
vorticity by integrating Eq. (39). 

The velocity field generated by the roll up process may be obtained once 
the location of the discrete concentrations of vorticity and the strength 
remaining in the unrolled sheets are known. Figure 8 illustrates this 
process. The velocity field at position "A" results from eight contributions: 
the right and left partially rolled up vortices, the right and left unrolled 
up sheets, and their image vortices and sheets. The velocity contribution 
from the right unrolled up sheet is computed from: 



( 41 ) 
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Figure 8. Schematic of how the velocity field at a particle position A is computed 




where: 
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y 2 + 


2 
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2 , . v 2 . 2 

r = (y - £) + z 


(42) 


and x is measured downstream of the wing and related to time by x = U^t • 
A similar expression gives the velocity contribution from the unrolled up left 
sheet and its image sheets. The average circulation is given by: 


r 

ave 



(43) 


where £ is the unrolled length. On the sheet the velocity is singular, 
since the sheet is infinitely thin. In reality the velocity varies linearly 
across the sheet. In the AGDISP code this linear variation is assumed to 
occur over a length equal to 0.l£ • 

The contributions of the rolled up portion of the vortices are similar. 
The velocity is of magnitude: 


U 


r(r) 
2tt r 


(44) 


where r is the distance from the centroid of the vortex to the particle 
position. The left and right vortices and their image vortices contribute to 
the magnitude of the velocity at the particle position by summing the 
contributions from Eq. (44). The velocity contribution of each vortex is 
perpendicular to a line connecting the vortex centroid and the particle 
position. 

In summary, the fixed wing wake flow field model is quite detailed and 
accounts for wake roll up and descent as well as for velocity field variations 
resulting from the details of the spanwise load distribution. The user may 
also elect to forgo this detail by specifying the position and strength of 
fully rolled up point vortices. 


Helicopter Flow Field 

The wake of a helicopter is considerably more complex than the wake of a 
fixed wing aircraft since the rotor downwash wake in hover is fluid 
dynamically different than the vortex wake generated by a fixed wing 
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aircraft. The flow field model included in the AGDISP code transitions from a 
hover wake at advance ratio y equal to zero, to a vortex wake at advance 
ratio equal to y Q . Figure 9 schematically shows the helicopter wake model 
option in the AGDISP code. The weight of the aircraft W is divided into a 
uniform downwash and a vortex wake according to: 



2p . ttR 2 U~ (h) 2 
air r 3 


(45) 


— = 2p . u r R 
y air 00 r r 


(46) 


where U^(h) is the induced velocity at the rotor computed from actuator disk 

theory and R r is the rotor radius. For values of y/y Q greater than unity, 

the complete weight of the helicopter is carried by the vortex wake. Eq. (46) 
defines the vortex strength F r of the tip vortices positioned at the end of 

the blade. At zero advance ratio the wake induced velocity at the rotor plane 
is uniform across the actuator disk, while at an advance ratio of y Q , the 

entire wake resembles that of a fixed wing aircraft. Beneath the rotor the 
downwash velocity, , from the actuator disk portion of the wake is assumed 

to decrease linearly to the surface as: 


(z) = U 3 (h)z/h 


(47) 


where h is the rotor altitude above the ground. 


Along the dividing streamline the pressure perturbation is zero, so that 
along this surface: 


2 2 2 

u 2 + u 3 = U 3 (h) 


(48) 


Since the flow is tangent to this surface: 
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dy s ~ U 2 


z 

s 


h 




(49) 


where the subscript s denotes the location of the dividing streamline. 
Integrating Eq. (49), subject to z g = h » Y s = R r yields: 





(50) 


for the y g position of the dividing streamline given the z g position of 

the dividing streamline. If a particle is located outboard of the dividing 
streamline, it will feel no effect of the helicopter downwash field. If a 
particle is within the downwash field, however, it will experience a downward 
velocity component given by Eq. (47) and a horizontal velocity component 
assumed to be linear from the wake centerline as: 


U 


2 



(51) 


When the advance ratio is zero the flow field is described by Eq. (47), (49) 
and (51) and the flow is contained in the initial plane shown in Figure 9 
(x = 0 since U,^ = 0). For non-zero advance ratios the initial plane moves 
downstream at velocity U^, and the agricultural material see the cross flow 
field resulting from a wake skewed downstream by the free stream. In 
addition, at non-zero advance ratios, a vortex pair flow field is added to the 
cross flow plane. The vortex strength T r is computed from Eq. (46) and the 
initial separation between vortices is the rotor diameter 2R r . 

The position of the centerline of the skewed wake is given by integrating: 


dz 


cr 

dx 



U 3 (h) 

U 


z 

h 


( 52 ) 


with z = h at x = 0 


The integration yields: 
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Propeller Swirl 

Fixed wing aircraft propulsion systems utilize a propeller driven by 
either a piston or turbo jet engine. The propeller flow field modeling in the 
AGDISP code treats the propeller as an actuator disk* At the disk the 
incremental velocity AU^ over the flight speed is related to the 

thrust T by: 


T = 2p irR AU, (U + AU.) 
air pl°° 1 


(54) 


where R is the propeller radius. In steady flight the thrust equals the 
drag so that: 


T = D 



P . u 2 s 
air 00 


(55) 


where S is the wing planforra area and C D is the aircraft drag 
coefficient. Combining Eqs. (54) and (55) to eliminate T , we obtain the 
induced velocity at the disk: 


AUj 1 + 

U = 




C D S 


(56) 


The propeller exerts a torque Q on the fluid resulting in an axial flux 
of angular momentum downstream of the actuator disk. If we assume that the 
distribution of swirling velocity is linear in distance from the propeller 
axis, then the axial flux of angular momentum may be computed as: 


Q = 



P . r 2 ft (U 
air p <= 


+ 


AU 1 )rdr 


T 0 a -:^rA^ U o O + Au .) 
Z air p p «> 1 


(57) 
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where ft is the angular velocity of the propeller. The useful power 
produced by a propeller is the thrust times the flight speed. This power is 
related to the actual power supplied to the hub through a propeller efficiency 
coefficient r\ : 


TU 

oo 


nQft 

q 


(58) 


where ft is the propeller rotational speed. Combining Eqs. (57) and (58) 
together with Eq. (55) yields: 


U 3 C S 
oo D 

P irnft R 4 (U + AU,) 
q p 00 1 


(59) 


showing that the swirl left in the propeller slip stream increases with 
decreasing propeller efficiency. 

Downstream of the actuator disk the propeller slip stream spreads 
primarily as a consequence of turbulent diffusion. At the actuator disk the 
mean squared turbulence level computed from superequilibrium theory is: 



( 60 ) 


where A is the turbulent integral scale and the gradient of slip stream 
velocity is estimated to be of the order of . If the largest eddies 

are assumed to be 60% of the radius of the slip stream (an assumption often 
used in jet flow calculations), the turbulence level at the actuator disk 
becomes: 


q 2 = 0.72AU 2 
P 1 


(61) 


By balancing the flux of turbulent kinetic energy in the slip stream with 
turbulent dissipation and assuming that the integral scale grows linearly with 
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downstream distance, the propeller wake area will grow as the square of the 
downstream distance. The turbulent kinetic energy downstream of the propeller 
in the slipstream can be shown to equal: 



(62) 


where x Q is the virtual origin computed from Ref. 20 as: 


x = 1.4U A/q 

o 00 p 


(63) 


Simple Terrain 

To investigate the influence of surface slope on deposition pattern, a 
simple terrain model is included in the AGDISP code. The interaction of the 
vortex wake with an inclined straight ground plane (Figure 10) is incorporated 
by modifying the position of the image vortex system to retain the requirement 
of zero fluid velocity perpendicular to the inclined plane. The helicopter 
and Betz roll up effects are similarly modified. 


Cross-wind 


Cross-wind flow can significantly alter the ground deposition pattern even 
when the aircraft is flying close to the surface. In a neutral atmospheric 
surface layer the horizontal velocity follows a logarithmic profile: 


U 2 (z) = ^(z^) in 



'ln(z Jz ) 
d o 


(64) 


where ^(z^) is the h° rizon tal velocity perpendicular to the flight path at 
a given altitude z. , and z Q is the surface roughness. This roughness is 
usually taken to be 1/30 of the actual physical roughness height.' 

Associated with this cross-wind profile is fluid turbulence resulting from 
the mean shear. With the assumption that the integral scale of turbulence is 
proportional to the distance from the surface (A = 0.65z) , the turbulent 
kinetic energy q becomes: 
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Figure 10. Vortex trajectories in the vicinity of inclined terrain 


(65) 
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The mean velocity logarithmic profile and its associated constant turbulence 
level are included in the AGDISP code. 


Plant Canopy 


The canopy flow field model programmed in the AGDISP code is an 
approximation of the second-order closure turbulence model of Wilson and 
Shaw^ Figure 11 illustrates the canopy configuration. The canopy is 

assumed to be of height h c , with a plant areal density denoted by A(z) . 
An effective canopy roughness height is defined as: 


h 

•/' 


A(z)zdz 



A(z)dz 


( 66 ) 


so that the mean cross-wind velocity above the canopy may again be logarithmic 
and given by: 


U 0 (z) - U 0 (z )£n(z/z )/ £n(z ,/ z ) , z > h 

Z Z d c d c c 


(67) 


The mean turbulence above the canopy is given by Eq. (65). The height z Q 
may be interpreted as the displacement thickness caused by the presence of the 
canopy. 

Within the canopy Itself the mean velocity and root mean square turbulence 
level are assumed to go to zero linearly as the ground (z=0) is approached. 
This approximation is consistent with the Wilson and Shaw data as seen in 
Figure 11. Thus for z < h c : 


U ? (z) =— U (z )£n(h /z )/£n(z /z ) 
z n Z d cc dc 


( 68 ) 
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Figure 11. Mean wind and root mean square turbulent fluctuations in a canopy (Wilson & Shaw 
Here U. = HL (z , ) /itn (z , /z ) 



q 2 = 0.845 
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The canopy not only modifies the turbulence and wind shear above and 
within the canopy but also results in a drag force on the fluid in the wake 
descending into the canopy. In the AGDISP code this drag force acts on any 
vortex entering the canopy (for simplicity, the helicopter wake resulting from 
the nonvortex portion of the downwash flow field is assumed to not interact 
with the canopy). Figure 12 illustrates the entrance of a vortex into the 
canopy. The retarding effect of the canopy is estimated by assuming that the 
swirling velocity in the entering vortex is reduced by a square drag law when 
entering the canopy: 




2 

2 

dt 


+ U„ 


(r,t) = - 


C c A(z)(U2 + 


u 2 ) 


(70) 


Since it is the circulation reduction that slows d own the vortex in the 
canopy, Eq. (70) is recast by substituting r = "bVu 2 + and defining a 
depth-averaged plant areal density as: 


A 



Adz 


(71) 


where a is the penetration depth of the vortex (as shown in Figure 12). If 
we assume that the retarding effect of the canopy will be proportional to that 
portion of the vortex flow field immersed in the canopy, Eq. (70) can be 
recast to become: 


d£ 

r 2 


C c AA'dt 

* 5 ^ 


(72) 


where A’/A^ is the fraction of the vortex within the canopy. At t = 0 we 
set r = r Q so that Eq. (72) yields: 
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Figure 12. Schematic of a vortex wake entering into a plant canopy 
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where 


t 



o 


( 74 ) 


Since F(t) is a monotonically increasing function of time, the effect of the 
canopy is to continually reduce the circulation of the vortex as it interacts 
with the canopy. The result of this reduction in circulation is a decrease in 
the descent rate of the vortex pair as it moves into the canopy, as well as a 
reduction in the swirling velocities and the wake convection of agricultural 
material. Consequently, a particle entering the canopy in the presence of 
cross-wind will feel a gradually reduced local velocity field tending to 
decelerate the particle. Vortices entering the canopy will become less 
effective in influencing particle motion because of the reduction in their 
circulation strength as given by Eq. (73). 


Superequilibrium Turbulence 

The option of including the detailed effects of turbulence on the 
dispersion of agricultural material is incorporated in the AGDISP code by 
invoking superequilibrium turbulent transport theory^ 22 '. Superequilibrium 
refers to the second-order closure turbulent transport model limit where the 
velocity correlations are able to track their equilibrium values. This limit 
is also referred to as the eddy viscosity limit, since velocity correlations 
are proportional to mean velocity gradients. For an incompressible flow the 
diffusion-free second-order correlations <u^Uj> satisfy: 
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Dissipation 



2 

where q = <u^u^> + <U 2 U 2 > = and repeated indices are summed. If the 

velocity correlations track the mean flow, D<u^Uj>/Dt = 0 , and: 
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6, . 3 

_Jj i. 

12 A 


= 0 
(76) 


Since the local mean flow gradients are known, the system of equation] 
represented by Eq. (76) can be solved exactly for <UjUj> to determine q‘ 
for the evaluation of particle dispersion. 
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4. APPLICABILITY OF AGDISP FOR LONG TIME DEPOSITION STUDIES 


The AGDISP code computes the motion of clusters of particles released from 
a single nozzle by predicting the ensemble averaged mean trajectory and 
variance about this mean. Since the flow field through which these particles 
pass is in general inhomogeneous, that is, the ensemble averaged mean velocity 
field, U i , and hence the mean square turbulence field, q 2 , are in general 
functions of position, there are restrictions with regard to spatial scales 
over which code predictions are accurate. 

A spatial scale defined by: 


L 


TS 



(77) 


may be seen to be a characteristic length over which the turbulent flow field 
changes significantly. Since the dispersion of the released material is 
predicted by assuming that fluid fluctuations acting on the cluster of 
particles are those fluctuations existing at the mean position, a simulation 
time may eventually be reached where variance of the cluster of particles will 
be on the order of L Tg . At this point an error is introduced in the AGDISP 
prediction of variance. 

Code predictions of dispersion are then reliable when the particle 
variance <x^x^> is such that: 


<x i x l > < l ts 


(78) 


Code predictions of variance may still be useful past this point, but the user 
should be aware of the approximation being made to estimate the variance. 
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5. AGDISP VALIDATION AND SAMPLE CASES 


Validation of the AGDISP code was undertaken by NASA and is documented in 
Ref. 23. Validation was made by comparing predicted deposition patterns against full 
scale measurements of deposition behind a fixed wing aircraft. To date, the code has 
not been validated for aerial application behind helicopters. 

A limited number of sample calculations are presented here to illustrate 
some of the deposition problems which may be investigated with the code. Note 
that in all calculations presented here the particle density is taken to be 
that of water. 


Helicopter with Boom Extending Outboard of a Rotor 

It is known that mounting a spray boom significantly outboard of the rotor 
may not necessarily increase swath width. A sample calculation is undertaken 
for a hypothetical helicopter with a gross weight of 13800 N and a rotor 
diameter of 8m. The spray boom is 6m above the ground, and the helicopter 
advance ratio is 0.2. The particles released from the boom have a specific 
gravity of unity and are 100 Mm diameter. 

Results of the computation are presented in Figure 13 where trajectories 
are shown for times up to 50 seconds after helicopter passage. It is clear 
that material released outboard of the rotor do not feel the favorable 
downwash of the rotor and remain adrift. The outward and initially upward 
motion of these outboard particles result from the portion of the helicopter 
wake which is vortex like since the advance ratio is non-zero. It is clear 
that the AGDISP code can be utilized to investigate optimal boom location and 
positioning of nozzles for a specific spray mission. 


Material Release in a Neutral Turbulent Atmospheric Boundary Layer 

In the AGDISP code the neutral atmospheric boundary layer has an assumed 
mean velocity profile which is logarithmic with altitude, while the 
atmospheric turbulence intensity is independent of altitude. Two computations 
are undertaken to illustrate the cross-wind effect. In the first (shown in 
Figure 14) particles of diameter 1000, 100 and 10 pro with specific gravity 
equal to unity are released from rest at an altitude of 6.2m. The mean wind 
at 6m is taken to be 0.6 m/sec with the surface roughness of 0.03m. This 
results in an rms turbulent velocity fluctuation of 0.1 m/sec independent of 
height in the neutral layer. As may be seen, the mean trajectory and 
dispersion are a very strong function of particle diameter. The dashed lines 
about the solid mean is the standard deviation of the particle computed normal 
to the trajectory. In Figure 15 particles of 100 \im diameter are released 
from rest at an altitude of 6.2m. The wind velocity at the release altitude 
is taken to be 0.6m and 1.5 m/sec. As anticipated the downstream point of 
impact and dispersion increase with higher wind speed. The greater dispersion 
results from the rms turbulent level increasing from 0.1 m/sec with a wind of 
0.6 m/sec at 6m to 0.25 m/sec with a wind of 1.5 m/sec at 6m. 
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Figure 13* Aerial application from a helicopter with spray boom extended outboard of rotor 
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Figure 14. Mean particle trajectories (solid) and standard deviation in a 
direction perpendicular to the mean (dashed). Particles of 
diameter 10 > 100 and 1000 ym released from rest at an altitude 
of 6 m . Mean wind at 6 m was 0.6 m/sec. 
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Figure 15. 


Effect of cross-wind magnitude on particle trajectory 
and on ground deposition pattern. 100 in particles are 
released from rest at an altitude of 6 m. Mean particle 
trajectory is shown solid and one standard deviation 
normal to the trajectory is shown dashed. 


Fixed Wing Aircraft with Triangular and Rectangular Spanwlse Load Distribution 


It is known that the spanwise load distribution of an aircraft can 
significantly alter the ground deposition pattern of material released into 
the wake. To illustrate this effect two load distibutions are considered: a 
wing which is triangularly loaded and a wing which is rectangularly loaded. 
To put the predictions of ground deposition on a consistent basis, it is 
assumed that aircraft weight, wing span, flight speed and altitude are the 
same between wing loadings. In the calculations undertaken here, wake roll up 
is neglected and it is assumed that the trailing vortices at the wing are 
located at the centroid of vorticity and downstream descend under their mutual 
influence. Referring to Figure 16 it may be shown** 1 ' that the swirling 
velocity of fully rolled up vortices from a triangularly loaded wing is 


u = 


u = 


2T 

T 

irb 

_It 

2tt r 


r < b/4 


r > b/4 


and 


(79) 


and from a rectangularly loaded wing is: 


1 

U « tt— for all r (80) 

2ttt 


where r is measured from y = b/2 at the wing elevation. In these 
calculations the wing span is 6.2m and particles are released from wing 
elevation at lateral positions which are 0, ±17, ±33, ±50, ±66 and ±83 percent 
of the semispan. It was assumed for convenience that the initial variance in 
particle position <x^x^> and velocity ^ v i v i^ (no SUTn ) were zero at the 

time of particle release. The required variances for fluid fluctuations were 
computed from the superequilibrium model described in Section 3. 

Results of the simulation with AGDISP are shown in Figures 17 and 18. The 
mean trajectory is shown as a solid line, while the dotted curves denote the 
magnitude of the variance of particle position computed normal to the 
trajectory. The effect of load distribution on 100 pm diameter particles is 
clearly seen. The higher wing root circulation in the triangularly loaded 
wing required to keep lift constant between computations results in particle 
trajectories being dominated by the vortices. In the rectangularly loaded 
case, only the particles released far outboard are dominated by the 
vortices. The results of these computations are, of course, not surprising 
since for a given particle size any reduction in the amount of swirling 
velocity in the wake is expected to diminish the influence of the wake on the 
deposition pattern. 
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Figure 16. Geometric details of the wake flow field 
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The particle trajectories and ground deposition pattern from a 
triangularly loaded wing after 10 seconds of simulation. A 
comparison of ground deposition pattern from a rectangularly 
loaded wing shown dashed is made with the deposition pattern 
from the triangularly loaded wing shown solid 




Cross-Wind 


It is well known that the effect of cross-wind has a significant impact 
with regard to deposition pattern and resulting drift. 

A calculation was undertaken with the horizontal velocity of 0.2 m/sec at 
6m altitude, with the assumed logarithmic profile shape. The surface rough- 
ness was taken to be 3 cm and the rms turbulence level was 10.7 cm/ sec. The 
particles were released at position ±y = ±0.38n m with n taking integer 
values between 0 and 13. The results of this computation is shown in Figure 
19 where the skewing of the deposition pattern on the ground is observed. The 
concentration distribution denoted by the dashed curve is the deposition 
pattern which results in the absence of cross-wind. It is noted that there 
results an excess in deposition upwind resulting from vortex swirling flow 
field countering the cross-wind. 


Evaporation 

The effect of evaporation is to reduce particle diameter and hence make 
particle dynamics more sensitive to wake and atmospheric flow fields. A 
computation was made with the wet-bulb depression A0 = 5 deg C identical to 
that illustrated in Figure 17 except that evaporation is permitted to occur. 
The results of the deposition computations are summarized in Figure 20 where 
it may be seen that when evaporation occurs there is nearly a factor of two 
reduction in mass deposited on the ground after 10 seconds. 
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Figure 19. Effect of cross-wind on particle trajectories and ground 

deposition pattern. Solid and dashed curves on the concen 
tration distribution are with and without a cross-wind 
respectively. The aircraft has a rectangular load distri- 
bution shown in Fig. 17. 







6. AGLINE 


The AGDISP code tracks the dynamics of material released from discrete 
locations. After diffusion has had sufficient time to act, a spray cloud is 
formed and material can no longer be identified with the location from which 
it was released. To incorporate this smearing feature of turbulence 
diffusion, an option in the AGDISP code, "AGLINE" ( Agricultural LINE 
Dispersal), can be exercised at any specified simulation time. This option 
replaces discrete distributions of material in the air with a single Gaussian 
distribution having the same mean and standard deviation as the discrete 
distributions. In this manner the AGDISP code can then be used to track the 
dynamics of a spray cloud or the AGLINE outputs can be used to initialize 
other dispersal codes. 

Use of the AGLINE feature is documented in the AGDISP User Manual. ^^) 
Figure 21 shows the results of using the AGLINE feature on the output of the 
sample calculation described in the User Manual. As can be seen in this 
example, four discrete distributions (dashed) are replaced by one distribution 
(solid) having the same mean and standard deviation. 
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Figure 21. Equivalent Gaussian distribution for the example test case in Ref. 16 




7. CONCLUSIONS AND RECOMMENDATIONS 


A computer code AGDISP has been developed which has the potential to 
reliably and efficiently predict the deposition of agricultural material 
released from rotary and fixed wing aircraft* This code computes the ensemble 
averaged mean motion of the material and the dispersion about this mean motion 
resulting from turbulent fluid fluctuations. These fluctuations result from 
turbulence generated by the aircraft itself or present normally in the 
atmosphere. Initial validation of the AGDISP predictions against deposition 
data taken behind a fixed wing aircraft is promising. Additional studies 
should be undertaken to confirm adequacy of the flow field modeling, 
particularly the helicopter wake downwash model, to assure acceptable accuracy 
of deposition predictions behind these aircraft. 
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